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ABSTRACT 
A  method  is  developed  to  analyze  and  compare  the  effectiveness  of 
ASW  vehicles.   The  measure  of  effectiveness  is  the  probability  that  the 
vehicle,  after  detecting  a  submarine  with  passive  sensors,  can  transit 
to  the  contact  area  and  re-establish  contact  with  the  submarine.   A  com- 
puter simulation  is  developed  and  an  example  using  three  hypothetical 
ASW  vehicles  is  illustrated. 
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1.    Introduction. 

As  a  result  of  increased  emphasis  currently  being  placed  on  ASW 
weapon  systems,  it  is  advantageous  to  possess  a  computer  simulation 
capable  of  analyzing  and  comparing  proposed  ASW  vehicles  and  opera- 
tional ASW  vehicles.   The  computer  simulation  for  the  comparison  of  ASW 
vehicles  developed  in  this  paper  is  intended  to  be  used  for  this  pur- 
pose. 

The  problem  at  hand  is  to  develop  a  method  for  comparing  ASW  vehi- 
cles under  operational  conditions  using  vehicle,  submarine,  and  tactical 
situation  parameters.   The  simulation  allows  one  to  compute  a  proba- 
bility which  is  used  as  a  measure  of  effectiveness.   This  probability, 
hereafter  referred  to  as  the  "probability  of  success,"  is  defined  in 
this  paper  as  the  probability  that  an  ASW  vehicle  successfully  relocates 
a  submarine  at  an  estimated  position  given  that  a  detection  has  occurred. 

A  description  of  the  simulation  and  its  limitations  is  made  in  Sec- 
tion 2.   It  is  very  important  to  understand  the  simulation  and  the 
background  assumptions  before  the  results  of  the  simulation  are  ana- 
lyzed.  After  the  problem  is  developed  mathematically  in  Section  2,  a 
discussion  and  an  illustration  is  given  in  Section  3  to  demonstrate  how 
to  use  the  computer  program.   In  Section  4,  several  applications  are 
discussed.   Finally  a  comparison  is  made  of  three  hypothetical  ASW 
vehicles--a  High  Speed  Surface  Vehicle,  a  VERTOL  Vehicle,  and  a  Seaplane- 
Type  Vehicle- -using  the  results  generated  from  the  computer  simulation. 
The  output  of  one  computer  run  is  illustrated  in  Figures  5  and  6.  Curves 
for  the  probability  of  success  versus  range  to  the  submarine  for  each 
vehicle  are  displayed  in  Figures  7,  8,  and  9.   These  curves  are  used  to 
analyze  and  compare  the  systems. 


The  results  of  these  comparisons  illustrate  the  usefulness  of  the 
simulation  in  comparing  the  performance  of  different  vehicles  and  for 
studying  the  effect  of  various  parameters  such  as  expected  sensor  range 
on  the  probability  of  success. 


2.   Simulation  Description  and  Mathematical  Development. 

Effort  has  been  directed  toward  the  development  of  a  computer  simu- 
lation that  represents  a  real  world  conflict  between  an  ASW  vehicle  and 
a  submarine.   The  setting  is  one  of  an  ASW  vehicle  assigned  to  a  patrol, 
surveillance,  or  barrier  mission.   Initial  detection  is  assumed  to  have 
been  made  using  passive  sensors.   The  vehicle  then  proceeds  as  rapidly 
as  possible  to  an  estimated  position  (EP)  in  order  to  relocate  the 
submarine.   The  EP  is  the  predicted  position  of  the  submarine  when  the 
actual  position  is  not  known,  but  the  vehicle's  sensors  indicate  the 
presence  of  a  submarine.   The  objective  of  the  vehicle  is  to  compute  an 
EP  of  the  submarine,  close  the  position  at  maximum  speed,  and  attempt 
to  relocate  the  submarine  using  active  sensors.   No  information  about 
the  actions  of  the  submarine  is  available  while  the  vehicle  closes  to 
the  EP. 

Certain  other  restrictions  are  assumed  in  the  model.   Ahead  thrown 
weapons  and  nuclear  weapons  are  not  available  for  use  by  the  vehicle. 
Multiple  targets  are  excluded.   A  passive  detection  range  of  at  least 
five  nautical  miles  is  necessary;  otherwise,  an  active  sensor  would  be 
used  and  the  problem  of  transiting  to  an  estimated  position  is  irrele- 
vant.  Also,  the  vehicle  is  assumed  to  lie  motionless  in  the  water  until 
detection  has  been  achieved  to  insure  passive  sensor  capability. 

The  simulation  is  developed  for  two  different  tactical  situations. 
The  first  case  is  that  of  a  submarine  which  is  unaware  of  the  activity 
of  the  vehicle.   The  second  case  is  that  of  a  submarine  which  is  alerted 
by  the  activity  of  the  vehicle.   Common  to  both  of  these  cases  is  a 
maneuvering  board  problem  which  will  be  discussed  first. 


The  maneuvering  board  solution  consists  of  describing  geometrically 
the  positions  and  motions  of  a  vehicle  and  a  submarine  with  time.   The 
passive  sensor  of  the  vehicle  provides  the  last  known  position,  course 
and  speed  of  the  submarine,  and  the  bearing  and  range  to  the  submarine 
from  the  vehicle.   Knowing  these  values  and  the  vehicle  speed,  an  EP, 
transit  time  of  the  vehicle,  relative  closing  speed,  and  the  distance 
traveled  to  EP  are  computed.   For  overall  ease  in  reading  the  same  nota- 
tion is  used  in  the  mathematical  development  that  has  been  used  in  the 
computer  simulation  (Appendix  C). 

The  solution  is  computed  by  solving  a  velocity  vector  problem.   The 
vehicle  is  located  at  the  origin  (0,0)  of  a  rectangular  coordinate  sys- 
tem.  One  side  of  the  velocity  triangle  (Figure  1)  represents  the 
estimated  velocity  vector  of  the  submarine  (Ul)  which  is  re-established 
at  the  origin.   The  velocity  vector  for  the  vehicle  (Z)  is  unknown  in 
direction,  but  known  in  magnitude.   The  remaining  relative  velocity  vec- 
tor (RELSPD)  represents  the  relative  course  and  speed  of  the  vehicle 
closing  the  submarine.   Only  the  direction  of  RELSPD  is  known  which  is 
parallel  to  the  true  bearing  line  from  the  vehicle  at  (0,0)  to  the  sub- 
marine at  (X2,  Y2).   Then  the  intersection  (X3,  Y3)  between  the  relative 
velocity  vector  and  the  vehicle  velocity  circle  is  computed. 

In  order  to  solve  the  maneuvering  board  problem  various  points  and 
vectors  in  the  rectangular  coordinate  system  must  be  computed.   The 
values  for  the  components  of  ul  are: 

XI  =  Ul  x  SIN  (THETA)  (1) 

Yl  =  Ul  x  COS  (THETA)  %  (2) 

The  last  known  position  (X2,Y2)  of  the  submarine  is  given  by: 
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X2  r  RAI  x  SIN  (BETA)  (3) 

Y2  =  RAI  x  COS  (BETA)  (4) 

where  BETA  is  the  bearing  in  degrees  of  the  submarine  from  the  vehicle 
and  RAI  is  the  range  in  nautical  miles  of  the  submarine  from  the  vehicle. 
The  RELSPD  vector  is  represented  by  the  line  through  the  point  (Xl,  Yl) 
and  parallel  to  the  line  connecting  the  origin  and  (X2,  Y2).   The  equa- 
tion of  this  line  is: 

Y3  =  ||  (X3-X1)  -+-  Yl  .  (5) 

The  RELSPD  vector  intersects  the  vehicle  velocity  circle  at  the  point 
(X3,  Y3).   The  equation  of  the  vehicle  velocity  circle  is: 

Y32  -+-  X32  =  V2^  (6) 

where  V  is  the  average  vehicle  transit  speed.   The  point  (X3,  Y3)  is  de- 
termined by  solving  equations  (5)  and  (6)  simultaneously.   The  direction 
of  the  vector  from  the  origin  to  (X3,  Y3)  represents  the  course  of  the 
vehicle  when  traveling  inbound  to  EP.  The  relative  speed  vector  is  rep- 
resented by  the  vector  between  (Xl,  Yl)  and  (X3,  Y3)  and  its  magnitude  is 

RELSPD  =  -W  (X3-Xl)2+  (Y3-Y1)2  .  (7) 

The  vehicle's  transit  time  (TD)  to  EP  is  defined  as: 

TD  =  RAI/RELSPD  .  (8) 

Let  the  vehicle  reaction  time  (RT)  be  the'  following  sum: 

RT  s  PTOT  +  TOT  +  CLOT  -I-  RELT  (9) 

where  PTOT  -  pretake-of f  time  which  commences  when  the  vehicle  can 
no  longer  observe  the  submarine  and  includes  decision 
time,  warm-up  time,  and  sensor  retrieval  time, 

TOT  -  take-off  time  which  includes  take-off  roll  and  time  to 
take-off  into  the  wind, 


CLOT  -  land  time  which  includes  time  to  land  and  taxi  time 

needed  to  arrive  at  EP  in  order  to  deploy  active  sensors. 

RELT  -  target  relocation  time  which  includes  the  time  necessary 
to  deploy  sensors  and  ends  when  the  submarine  is  re- 
located. 

Any  of  these  times  can  have  a  zero  value.   For  example,  a  helicopter 

could  have  TOT  equal  to  zero  if  it  needs  no  take-off  time.   Also,  the 

times  are  mutually  exclusive  so  that  no  time  interval  is  added  twice. 

Obviously  the  blind  time  (TB)  is: 

TB  '   RT  -f-  TD.  (10) 

Evasion  for  the  alerted  submarine  begins  when  blind  time  commences.   At 
that  time  the  submarine  is  assumed  to  be  at  (X2,  Y2).   If  RT  is  equal  to 
zero,  the  submarine  is  at  (X2,  Y2)  when  evasion  begins;  but  if  RT  is 
not  equal  to  zero,  the  submarine  is  not  actually  at  (X2,  Y2)  when  eva- 
sion begins  because  it  takes  RT  for  the  vehicle  to  commence  moving 
towards  EP.   During  the  reaction  time  the  vehicle  is  motionless,  but 
the  submarine  is  moving  away  from  (X2,  Y2).   Because  of  the  small  dis- 
placement of  the  submarine  from  (X2,  Y2)  due  to  RT,  the  assumption  that 
the  submarine  is  at  (X2,  Y2)  when  evasion  begins  does  not  significantly 
affect  the  results  of  the  simulation. 

It  is  necessary  to  determine  the  estimated  position  (EP)  of  the 
submarine  after  a  time  interval  of  TB  and  a  submarine  speed  U3.  The 
equation  of  the  estimated  distance  (DIST)  traveled  by  the  submarine  is: 

DIST  =  U3  x  TB  .  (11) 

The  EP  is  also  the  position  at  which  the  vehicle  estimates  that  an  in- 
tercept will  be  made  with  the  submarine.   The  components  of  EP  relative 
to  the  origin  are: 

X4  =  DIST  x  SIN  (THETA)  +  X2  (12) 


Y4  =  DIST  x  COS  (THETA)  -+-  Y2  .  (13) 

The  equation  for  the  distance  DISDAT)  traveled  by  the  vehicle  to  EP  is 

DISDAT  =  ^  X42  +  Y42  .  (14) 

The  values  for  TB  and  DIST  are  very  important  in  determining  the  proba- 
bilities of  success  for  the  two  cases  developed  below. 

Case  I.  The  unalerted  submarine: 

It  is  conceivable  that  the  vehicle  could  transit  to  EP  without  the 
submarine  being  cognizant  of  the  impending  danger.   In  this  case  the  un- 
alerted submarine  does  not  evade  and  proceeds  at  cruise  course  and  veloc- 
ity.  Prediction  errors  for  the  submarine  course  and  velocity  are  made 
which  are  characteristic  of  the  vehicle's  equipment.  The  errors  are 
assumed  to  be  distributed  normally  about  the  predicted  course  and  veloc- 
ity.  Because  of  long  blind  times  expected,  the  predicted  positions  of 
the  submarine  are  not  accurately  represented  by  a  circular-normal  distri- 
bution, and  a  numerical  integration  must  be  performed.  Probabilities  of 
success  are  computed  for  each  of  a  representative  distribution  of  pre- 
dicted positions  consisting  of  25  points  or  cells  (see  Figure  2).   This 
concept  has  been  adopted  from  \3 J . 

The  last  known  submarine  position  (X2,  Y2)  is  established  at  the 
origin  of  a  rectangular  coordinate  system.   Since  DIST  is  the  distance 
traveled  by  the  submarine  during  TB,  then  EP  is  (DIST,0)  in  this  coor- 
dinate system. 

Using  the  prediction  errors,  a  submarine  position  is  computed  for 
each  of  the  25  cell  midpoints.   The  distance  (W)  and  angle  (THETER)  for 
each  cell  are  computed  using  the  submarine  cruise  speed  (Ul),  the  submarine 


8 


SBIAS(I,T) 


O^UERR 
<r2=THDEL 


Figure.  2.  Numerical     I«Vea*"eu"»ot\     For  TV.t.     U«aler"fa4     C*s«» 


speed  error  (UERR),  the  submarine  course  error  (THDEL) ,  and  the  blind 

time  (TB).   The  equations  for  W  and  THETER  are: 

W  -    (Ul-AA(I)  x  UERR)  TB  (15) 

THETER  =  AA  (J)  x  THDEL  ,  (16) 

with  AA(I)  and  AA(J)  *  0,  +  .52,  ±  1.28  depending  on  which  cell  is  being 

considered.   These  numbers  can  be  obtained  from  any  normal  density  tables, 

K. 

The  equations  for  the  submarine  position  (SUBX  (I, J),  SUBY(I,J))  at 

the  center  of  each  of  the  25  cells  are: 

SUBX  (I, J)  «  VV  x  COS  (THETER)  (17) 

SUBY  (I, J)  =  W  x  SIN  (THETER)  .  (18) 

A  bias  is  computed  between  EP  and  the  center  of  each  cell.  The  equation 

for  the  bias  of  the  (i,j)th  cell  is: 


SBIAS  (I, J)  =  -J  (SUBX(I,J)-DIST)2+-  (SUBY(I,J))2  .         (19) 
Various  errors  have  to  be  considered  to  give  the  simulation  realism. 
It  is  assumed  that  the  navigation  error  and  sensor  locating  error  are 
additive  about  the  point  EP.   The  distribution  of  the  vehicle  location 
about  EP  due  to  these  errors  is  assumed  to  be  circular  normal. 
The  standard  deviation  of  navigation  error  (SIGNAV)  is: 

SIGNAV  =  PRCNTN  x  RAI  (20) 

where  PRCNTN  is  a  percentage  of  the  last  known  range  (RAI).   The  stand- 
ard deviation  of  sensor  locating  error  consists  of  the  standard  devia- 
tion of  bearing  (  SIGSBR  )  and  the  standard  deviation  of  range  (SIGSR) 
which  are  defined  by: 

SIGSR  =  PRCNTS  x  RAI  (21) 

SIGSBR  *  RAI  x  SIN  (BRGER)  (22) 

10 


where  BRGER  is  the  bearing  error  and  PRCNTS  is  the  percentage  error  of 
the  last  known  range.   Since  (21)  is  greater  than  (22),  because  of  the 
passive  sensor,  an  elliptical  normal  distribution  exists;  but  it  is 
approximated  by  a  circular  normal  distribution  for  ease  of  computation. 
This  is  accomplished  by  equating  the  area  of  an  ellipse  to  the  area  of  « 
circle.  The  result  is  an  approximate  standard  deviation 


XE    YE         CIR        u  XE  =  SIGSBR  (23) 

0"ye  =  SIGSR 

.'.    SIGAPR  =  #"CIR  ="yO"xE  0~YE  (24) 

Therefore,  the  total  standard  deviation  of  error  (SIGJT)  is  the  sum 
of  the  squares  of  the  approximate  standard  deviation  of  error  and  the 
standard  deviation  of  navigation  error  according  to  the  equation: 


SIGJT  *ll  (SIGNAV)2-H  (SIGAPR)2   .  (25) 

Other  error  models  were  considered  in  order  to  more  accurately  represent 
the  errors  involved.  Due  to  lack  of  time  the  error  model  described 
above,  although  not  the  best,  was  adopted. 

To  compute  the  probability  of  success  for  the  unalerted  submarine 
the  circular  coverage  function  is  used, \l\ .   In  order  to  use  this  func- 
tion, the  parameters  SBIAS(I,J),  SIGJT,  and  EXPRNG  are  required.   The 
first  two  have  been  defined  above.   The  term  EXPRNG,  which  is  defined  as 
the  expected  range  of  the  active  sensor  used  by  the  vehicle,  is  a  charac- 
teristic of  the  sensor  employed.   The  quantities  SBIAS(I,J)  and 

SIGJT 

EXPRNC 

-  are  used  in  the  circular  coverage  function  to  compute  the  proba- 
SIGJT 

bility  of  success  (PK(I,J))  that  a  circular  disk  of  radius  EXPRNG  will 
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cover  a  point  (SUBX(I.J),  SUBY(I,J)  from  the  EP  since  the  probable  posi- 
tion of  the  disk  is  described  by  a  Gaussian  distribution.   This  proba- 
bility PK(I,J)  is  the  probability  of  success  for  the  ijth  cell.   Since 
it  is  equally  likely  that  the  submarine  is  in  each  of  the  25  cells, 
the  probability  that  it  is  within  each  cell  is  .04.   The  probabilities 
of  success  for  each  cell  are  summed  and  multiplied  by  .04  to  give  the 
"probability  of  success  (PROB)  against  an  unalerted  submarine."  The 
probability  of  success  is: 


PROB  =  <^   >     (PK(I,J)  x  .04)  (26) 

J 

Case  II.   The  alerted  submarine: 

The  submarine  is  aware  of  the  presence  of  an  ASW  vehicle  and  evades 
by  changing  course  and  speed.   The  speed  changes  are  characteristics  of 
the  submarine  but  the  course  changes  are  determined  in  number  and  degree 
by  the  situation  being  simulated.   A  probability  of  success  is  computed 
for  each  evasion  turn  and  an  arithmetic  average  taken  over  the  number 
of  evasion  tactics  used. 

This  case  is  designed  to  give  practically  full  evasion  to  the  sub- 
marine (see  Figure  3).   The  submarine  evasion  turns  are  limited  to  a 
maximum  turn  of  90  degrees  to  either  side  of  the  estimated  submarine 
course.   No  additional  turns  are  granted  after  the  initial  turn  is  exe- 
cuted.  For  turns  greater  than  90  degrees  the  submarine  would  deviate 
from  its  initial  track  to  such  an  extent  that  it  would  not  be  able  to 
make  up  the  lost  time  necessary  to  make  good  a  mission  speed  of  advance. 

The  range  of  values  for  the  set  of  evasion  turns  is  determined  by 
the  situation  being  simulated.   For  example,  a  slow  submarine  would 
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almost  always  turn  away  from  the  vehicle  when  the  vehicle's  presence  is 

known.   There  are  two  limiting  turns  which  bound  the  set  of  evasion  turns 

to  the  right  and  the  left.   The  bounds  must  be  no  greater  than  90  degrees 

to  either  side  of  the  estimated  submarine  course.   PSI  and  DELTA  are  the 

left  bound  and  right  bound,  respectively.   The  range  of  values  for  the 

set  of  evasion  turns  is  determined  by  (PSI-DELTA).   The  interval  (ETA) 

between  evasion  turns  is  ETA  =  PSI-DELTA  f   NU  being  the  number  of  eva- 

NU 
sion  turns  considered.   If  the  set  of  turns  encompasses  turns  to  the 

right  of  the  submarine  course  then  ETA  =  ?-§! — (360-DELTA)  §   In  Figure  3 

NU 

PSI  =  060  degrees,   DELTA  =  330°,  and  ETA  =  060  (360-330)  =  15°. 

6 

Therefore,  the  evasion  turns  considered  for  this  example  are  330°,  345°, 
360°,  015°,  030°,  045°,  and  060°. 

In  the  alerted  submarine  case  the  submarine  course  and  velocity  pre- 
diction errors  that  were  present  in  the  unalerted  case  are  ignored. 
Since  in  this  case  a  set  of  evasion  turns  are  selected  which  are  only  sub- 
jective predictions,  it  is  valid  to  exclude  prediction  errors  because 
they  are  comparatively  insignificant  to  the  computations.  However,  ve- 
hicle sensor  and  navigation  errors  that  were  assumed  in  the  unalerted 
case  still  prevail. 

Initially,  the  submarine  is  at  cruise  speed.   When  evasion  commences, 
the  assumed  submarine  speed  is  increased  to  fast  speed.   The  character- 
istics of  the  submarine  determine  how  long  the  submarine  remains  at  fast 
speed  (length  of  time  is  Tl)  after  which  it  slows  to  silent  speed.   The 
time  intervals  of  acceleration  and  deceleration  are  considered  negligible 
when  compared  with  the  long  blind  times  encountered. 

The  alerted  case  uses  a  rectangular  coordinate  system  with  the  last 
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known  position  of  the  submarine  (X2,Y2)  at  the  point  of  origin  when  eva- 
sion begins.   From  this  point  positions  (X(N),  Y(N)),  which  are  the 
positions  at  the  end  of  each  of  the  N  evasion  turns,  are  calculated. 
(See  Figure  3.)  The  time  to  turn  (TTT),  evasion  turn  in  degrees  (AB), 
and  the  turning  radius  (TR)  are  needed  to  compute  X(N)  and  Y(N).   The 
equations  for  TTT,  X(N) ,  and  Y(N)  are  computed  as  follows: 

.p™,  _  arc  length  _  TR  x  AB 
sub  speed      U3 

If  TI  ^  TB,  then 

X(N)  »  TR  x  SIN(AB)-MJ3  x  (TB-TTT)  COS(AB)  (27) 

Y(N)  =  TR-  TR  xC0S(AB)-hU3x  (TB-TTT)  SIN(AB)  .  (28) 

If  TI  <  TB,  then 

X(N)  =  TRx  SIN(AB)-hU3x  (TB-TTT)  COS(AB)+  (TB-Tl)U2xCQS(AB)  (29) 

Y(N)  »  TR~TRxCOS(AB)+  U3  x (TB-TTT)SIN(AB) + (TB-Tl) 

U2x  SIN(AB).   (30) 

The  time  Tl  is  the  length  of  time  the  submarine  remains  at  fast  speed  (U3); 

after  which  it  changes  to  silent  speed  (U2).   Equations  (29)  and  (30) 

account  for  the  speed  change  occurring  before  TB  has  terminated.   The 

point  (X(N),Y(N))  can  be  considered  the  submarine's  actual  position  after 

the  N   evasion  turn  is  taken. 

The  estimated  position  EP  is  the  predicted  position  of  the  submarine 

after  the  interval  of  blind  time  has  elapsed.   Therefore,  EP  is  (DIST,0). 

The  distance  between  each  (X(N),Y(N))  and  EP  constitutes  a  bias  (BIAS(N)). 

This  is  an  important  factor  in  computing  the  probability  of  success.   The 

bias  for  the  Ntn  evasion  turn  is: 


BIAS(N)  =  -J   (X(N)-DIST)2+-  (Y(N))2  .  (31) 

As  in  the  unalerted  case,  the  parameters  BIAS,  EXPRNG,  and  SIGJT  must  be 
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computed  (EXPRNG  and  SIGJT  are  the  same  as  in  the  unalerted  case)  in 
order  to  use  the  circular  coverage  function.   The  circular  coverage  func- 
tion is  used  to  compute  the  probability  (PPK(N)),  for  each  N,  that  the 
vehicle  can  regain  contact  at  EP  given  the  total  standard  deviation  of 
error  (SIGJT),  the  bias  (BIAS(N))  between  the  terminal  point  of  the  Nth 
evasion  and  EP,  and  the  expected  active  sonar  range  (EXPRNG). 

The  two  ratios  of  the  parameters  that  are  needed  for  the  circular 

coverage  function  are  BIAS(N)  and  EXPRNG  t  with  these  ratios  the  proba- 

SIGJT       SIGJT 

bilities,*  PPK(N),  are  computed  for  each  N,  and  then  all  PPK(N) ,  N  ■  1, 
...,  NU+1,  are  summed  and  averaged  over  all  NU  +  1.   The  average  proba- 
bility is: 

NU+1 


PPK(N) 
PAVG  =  tTTT    __  (32) 

NU  +  1 
It  might  be  appropriate  to  use  the  option  to  weight  each  of  the  eva- 
sion turns.  If  it  is  believed  that  the  turns  within  the  set  of  evasion 
turns  are  not  equally  likely,  then  the  probability  vector  UNU(N) ,  N"l, 
. . . ,  NU+  1  can  be  used  which  gives  probabilities  that  each  of  the  N 
evasions  is  actually  taken.   If  they  are  not  equally  likely,  then  the 
"weighted  probability  of  success  (WTPAVG)  against  an  alerted  submarine*1 
is  computed.   The  equation  for  WTPAVG  is 

NU  +  1. 

^T-    (PPK(N)  x  UNU(N)) 
WTPAVG  =  ^fTT  (33) 

(NU+-1) 

The  results  of  the  mathematical  development  have  produced  methods 
of  computing  the  probability  of  success  (PAVG)  for  an  alerted  submarine 


16 


with  equally  likely  evasion  turns,  the  probability  of  success  (WTPAVG) 
for  an  alerted  submarine  with  evasion  turns  that  are  not  equally  likely, 
and  finally  the  probability  of  success  (PROB)  for  an  unalerted  submarine, 
These  probabilities  are  used  to  analyze  and  compare  ASW  vehicles. 
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3.    Input  Data  Rules. 

A  complete  run  for  the  computer  program  (Appendix  C)  consists  of 
25  input  cards.   Each  input  data  card  represents  the  value  of  one  vari- 
able except  for  the  last  two  cards,  one  of  which  contains  the  values  for 
the  AA(i),  i-1,  ...,  5  and  the  other  which  contains  the  UNU(N)  values, 
N»l,  ...,  NU+1  . 

One  of  the  first  functions  of  the  program  is  to  set  the  indices 
N0=1  and  Msl.   (See  Appendix  B--Flow  Charts.)   The  index  NO  is  used  to 
number  the  pages  of  the  program  output.   The  index  M  is  compared  with 
NEXT  (initially  set  equal  to  one).   The  comparison  dictates  whether  the 
program  will  continue  or  terminate.   The  value  for  NEXT  determines  the 
number  of  runs  to  be  executed  and  is  an  input  to  the  program.   The  compu- 
ter then  reads  the  input  data  cards  from  SUBROUTINE  INPUT  after  which 
the  25  values  are  printed  on  page  1  of  the  output.   The  index  NO  is  in- 
creased by  one.   The  programmed  computations  are  performed  and  the 
results  are  printed  on  page  2  of  the  output  (Figure  6).   The  indices  NO 
and  M  are  then  increased  by  one.   The  index  M  is  compared  to  NEXT  and 
if  NEXT  is  greater  than  M,  SUBROUTINE  INPUT  searches  for  more  input  data 
cards.   Since  the  program  is  used  to  compare  and  evaluate  systems,  only 
a  few  variables  will  be  changed  for  each  successive  run.   However,  more 
or  all  the  variables  can  be  changed  as  long  as  the  input  procedure  is 
followed.   For  example,  to  compare  the  three  vehicles  in  Table  1,  first 
place  the  parameters  for  the  High  Speed  Surface  Vehicle,  RAI  =  5.0,  and 
a  value  for  NEXT  in  the  data  input  for  Run  1.   (See  Figure  4.)   The  next 
run  consists  of  a  Locator  Card  which  initiates  the  reading  of  Variable 
Change  Card,  which  in  turn  changes  the  value  of  RAI  to  10.   The  subsequent 
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card  would  be  the  Information  Card,  which  indicates  no  further  changes 
for  the  second  run  --  start  computations.   For  this  run  the  variable  RAI 
has  changed  and  all  other  parameter  values  remain  the  same. 

If  another  variable  change  were  desired,  another  Locator  Card, 
Variable  Change  Card,  and  Information  Card  would  be  placed  after  the 
RAI  Variable  Change  Card  as  indicated  in  Run  3  in  Figure  4. 

A  run  following  the  initial  run  takes  a  minimum  of  three  cards.   If 
only  one  variable  is  being  changed  (e.g.  RAI  =  10),  then  three  cards  are 
necessary.   This  procedure  is  followed  for  RAI  a   5,  10,  ...,  60  for  each 
system  in  Table  1  in  order  to  compute  probabilities  of  success  for  each 
vehicle  over  the  entire  RAI  range. 
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4.    An  Example  and  Some  Applications. 

The  simulation  possesses  the  capability  of  being  used  as  a  tool  to 
analyze  and  compare  ASW  vehicles  as  explained  in  the  previous  sections. 
In  order  to  illustrate  this  capability,  the  parameters  listed  for  each 
system  in  Table  1,  and  ranges  RAIa5,  10,  15,  ...,  60  nautical  miles 
were  used  as  input  data  in  the  computer  program  (Appendix  B)  and  pro- 
cessed in  the  CDC-1604  computer.   The  program  is  written  in  FORTRAN-60. 
The  flow  chart  is  illustrated  in  Appendix  C. 

A  sample  output  for  the  Seaplane  Type  Vehicle  using  the  parameters 
of  Table  1  and  RAI"45  is  shown  in  Figures  5  and  6.   The  former  is  a  list 
of  all  the  parameter  values  for  the  indicated  run.   This  list  enables 
one  to  insure  that  the  correct  parameters  for  that  run  are  being  properly 
read  into  the  computer.   The  latter  page  contains  the  probabilities  of 
success  for  one  run.   The  probabilities  for  each  RAI  are  plotted  on 
graphs.   Each  vehicle  is  depicted  on  an  individual  graph.   Each  graph 
has  six  curves  --  three  drawn  from  values  calculated  from  using  a  slow 
submarine  as  an  adversary  and  three  drawn  from  using  a  fast  submarine  as 
an  adversary  (see  Figures  7,  8,  9). 

The  parameter  values  in  this  sample  are  strictly  hypothetical  and 
demonstrate  only  the  use  of  the  simulation.   However,  with  these  chosen 
parameter  values  it  can  be  seen  in  Figures  7,  8,  and  9  that  throughout 
the  range  of  RAI  the  High  Speed  Surface  Vehicle  has  probabilities  of 
success  considerably  lower  than  the  other  systems.   The  slower  the  speed 
of  the  vehicle,  the  lower  the  probability  because  of  a  low  closing  rate 
and  long  blind  time.   The  probabilities  of  success  are  higher  at  shorter 
RAI  for  the  VERTOL  Vehicle  as  compared  with  the  Seaplane-Type  Vehicle, 
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but  a  crossover  point  occurs  at  approximately  RAI=25,  where  the  Seaplane- 
Type  Vehicle  has  higher  probabilities  of  success.  The  crossover  is  a 
result  of  the  values  selected  in  the  sample  for  the  vehicle  speed,  reac- 
tion time,  and  expected  sonar  range.   From  these  results  the  best  vehicle 
can  be  chosen  and  the  other  vehicles  investigated  for  improvements. 

In  addition  to  the  example  stated  above  several  applications  are 
suggested  for  the  simulation.   One  application  is  to  use  the  simulation 
to  determine  what  parameters  or  characteristics  the  vehicle  and  the 
sensors  must  have  to  obtain  a  given  probability  of  success.   For  example, 
if  the  state-of-the-art  for  active  sensors  dictates  that  five  nautical 
miles  is  the  maximum  range  expected  for  active  sensors  for  the  next  ten 
years,  then  it  would  be  advantageous  to  find  values  of  vehicle  speed, 
reaction  time,  navigation  errors,  and  prediction  errors  that  can  be  tol- 
erated for  the  vehicle  and  still  maintain  a  given  probability  of  success. 
The  development  of  a  new  ASW  vehicle  could  be  directed  towards  achieving 
these  parameter  values. 

Another  application  can  be  made  if  a  new  sensor  is  being  proposed 
for  a  given  vehicle.   The  values  for  the  expected  active  sensor  range 
and  the  sensor  errors  can  be  computed  for  a  given  probability  of  success. 
These  parameter  values  can  be  used  as  design  limits  in  the  development 
of  the  proposed  sensor. 

Finally  the  simulation  can  be  used  if  intelligence  indicated  that 
a  new  enemy  submarine  is  in  existence.   The  simulation  can  be  used  to 
indicate  how  well  the  present  vehicles  will  perform  against  the  new 
threat  relative  to  the  old  threat. 
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5.    Conclusions. 

The  example  in  Section  4  is  the  type  of  comparison  that  can  be  per- 
formed for  various  kinds  of  ASW  vehicles.   There  are  several  dec is ions -- 
evasion  tactics  used,  submarine  type  used,  and  so  on- -that  have  to  be 
made  prior  to  using  the  program;  but  once  the  decisions  are  made  for  a 
given  situation,  they  have  to  remain  the  same  for  all  systems  being  com- 
pared.  Several  modifications  and  suggestions  for  the  model  are  discussed 
below. 

It  might  be  argued  that  an  elliptical  coverage  function  should  be 
used  instead  of  a  circular  coverage  function  because  of  the  elliptical 
sensor  locating  error.   This  change  can  be  made  by  rewriting  SUBROUTINE 
OCIP  of  the  computer  program  (Appendix  C). 

Instead  of  evasion  turns  limited  to  90  degrees  either  side  of  the 
course  line,  they  could  be  programmed  for  360  degrees  evasion;  i.e.,  180 
degrees  either  side  of  the  course  line. 

It  must  be  realized  that  the  depth  of  the  submarine  has  no  direct 
influence  on  this  simulation.   Therefore,  this  has  to  be  kept  in  mind 
when  parameter  values  such  as  EXPRNG  are  put  into  the  simulation.   It  is 
possible  that  EXPRNG  could  be  less  in  an  evasion  problem  where  the  sub- 
marine is  more  likely  to  increase  its  depth,  and  therefore,  EXPRNG  should 
be  reduced. 

It  would  be  possible  to  give  the  vehicle  freedom  of  movement  before 
transit  commences,  rather  than  being  motionless.   But  the  advantages  of 
the  passive  tactic  for  the  vehicle  are  then  relinquished. 

The  restriction  on  Tl  in  the  alerted  submarine  case  could  be  abol- 
ished by  programming  for  the  possibility  of  having  a  speed  change  exe- 
cuted before  a  turn  is  completed. 
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Possibly  a  passive  sensor  detached  from  the  vehicle  could  transmit 
current  information  about  the  submarine's  actions  while  the  vehicle  is 
inbound  to  intercept  the  submarine.   This  would  increase  the  probability 
of  success. 

In  the  foregoing  sections,  the  simulation  has  been  described  and  the 
mathematical  development  shown.   In  addition,  the  use  of  the  computer 
program  has  been  illustrated  and  an  example  and  its  results  briefly  ana- 
lyzed.  The  use  of  the  simulation  is  simple  and  the  results  easily  under- 
stood.  Therefore,  it  is  recommended  that  the  use  of  this  simulation  be 
employed  to  compare  and  analyze  proposed  ASW  vehicles  and  operational 
ASW  vehicles. 
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Notation 
A(1),...,A(5) 


APPENDIX  A 

GLOSSARY 

Meaning 

Number  of  standard  deviations  from  the  mean  of 
the  normal  density  function.  Each  A(i)  is  the 
midpoint  of  a  segment  of  area  which  represents 
20%  of  the  total  area  under  the  normal  density 
function.  The  A(i)'s  are  used  in  the  predic- 
tion error  of  the  unalerted  submarine  case. 


AB 


Evasion  turn  in  degrees  from  the  estimated 
submarine  course. 


BETA 


BIAS(N) 


True  bearing  of  the  submarine  from  the  vehicle 
(degrees) . 

Distance  from  EP  to  the  terminal  point  of  the 
jjth  evasion. 


BRGER 
CLDT 

DELTA 

DIST 

EP 

EXPRNG 

NEXT 
NU 

PAVG 


Sensor  bearing  error  in  degrees. 

Land  time  (min.)  which  includes  time  to  land 
into  the  wind  and  taxi  time  to  get  into  posi- 
tion to  employ  active  sensors. 

Bound  to  the  right  of  the  known  submarine 
course.  This  bounds  to  the  right  the  set  of 
possible  evasion  turns  (degrees). 

Distance  the  submarine  travels  during  blind 
time,  assuming  a  constant  course  and  speed. 

Estimate  position  of  the  submarine  after  blind 
time  has  expired. 

Expected  range  of  the  active  sensor  used  by  the 
vehicle  in  re-establishing  contact  at  the  EP 
(in  nautical  miles). 

-  Number  of  data  runs  for  the  computer  program. 

Number  of  increments  considered  with  the 
(PSI-DELTA)  range  of  the  set  of  possible  turns. 

Average  probability  of  success  for  N  evasions 
considered  in  the  alerted  submarine  case. 


PRCNTN 


Navigation  error  in  percent  of  distance  trav- 
eled. 
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Notation 

PRCNTS 

PROB 

PSI 
PTOT 


RAI 
RELT 

RT 

SBIAS(I,J) 
SIGJT 


Meaning 

Sensor  range  error  in  percent  of  sensor  range. 

Probability  of  success  for  the  unalerted  sub- 
marine case. 

Bound  to  the  left  of  the  last  known  submarine 
course.   This  bounds  to  the  left  the  set  of  pos- 
sible evasion  turns  (degrees). 

Pre  take-off  time  (min.)  which  commences  when 
the  vehicle  can  no  longer  observe  the  submarine 
and  includes  decision  time,  warm-up  time,  and 
sensor  retrieval  time. 

Range  from  the  vehicle  to  the  submarine  at  the 
last  known  position  (nautical  miles). 

Target  relocation  time  (min)  which  includes  the 
time  to  deploy  sensors  and  ends  when  the  sub- 
marine is  relocated. 

Reaction  time  which  is  the  sura  of  PTOT,  TOT, 
CLDT,  and  RELT. 

Distance  between  the  (i,j)   cell  and  EP. 

Total  standard  deviation  of  errors. 


TB 


TD 


Tl 


THDEL 
THETA 

TOT 

TR 


Blind  time  which  is  the  sum  of  the  reaction 
time  (RT)  and  the  transit  time  (TD). 

-  Time  it  takes  the  vehicle  to  travel  from  its 
initial  position  to  EP. 

-  Tactical  interval  of  time  (min.)  after  evasion 
commences  during  which  the  submarine  travels  at 
fast  speed.   After  Tl  has  elapsed,  the  sub- 
marine slows  to  silent  speed. 

Submarine  course  error  (degrees). 

-  Course  of  the  submarine  (degrees)  at  the  last 
known  position. 

-  Take-off  time  (min.)  which  includes  take-off 
roll  and  the  time  to  maneuver  into  the  wind. 

-  Turning  radius  of  the  submarine. 
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Notation  Meaning 

UERR  -   Submarine  speed  error  (knots). 

UNU(1 ),...,  UNU(NU+1)  -  The  set  of  probabilities  for  N=l,...,  NU+1 

evasion  turns.   These  are  used  when  evasion 
turns  considered  are  not  equally  likely. 

Ul  -  Submarine  cruise  speed  (knots). 

U2  -  Submarine  silent  speed  (knots). 

U3  -  Submarine  fast  speed  (knots). 

V  -  Average  vehicle  transit  speed  (knots) 

WTPAVG  -  Weighted  probability  of  success  for  the  alerted 

submarine  case. 

(X2,Y2)  -  Last  known  position  of  the  submarine. 

(X(N),Y(N))  -  Terminal  point  of  the  Nth  evasion  after  blind 

time  has  expired. 

Z  -  Average  vehicle  transit  speed  (knots).   Neglect 

acceleration  and  deceleration  time. 
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APPENDIX   B 


NIO  =  l 


3 


NEYT-M 


\T>o 


CALL 


INPUT 


<o 


INPUT  AND 
UNIT  CHANGES 


W   STOP j 


INPUT 
PARAMETERS 


CHANGE  SPEEDS 

FROM    KNOTS 

TO 

NAUTICAL  MILES 

PER,  MINUTE 


CHANGE  ANGLES 
FROM  DEGREES 
TO 
RADIANS 


3' 5 


MANEUVERING  BOARD 
PROBLEM 


Y>UU$IN(THETA) 

Yl=Ul  *cos(theta) 


X"2=RAIxSiw(G£TA.) 
Y2.»RMxCOS(6£T/\> 


GEDi 


X3(0=f\MYl)1 
X3fc)»-K3(0 


© 


ABlXX=|xi-X3(>)| 
AB2XX=|X^-X3(Z)| 


=  0 


ASIXX-ASIXX 


>0 


\/_ 


ERROR 
STATEMENT 


<0 


V3(I)=Y| 


Xi'O 


YXO=Vvz-(x')2 
Y3ft)=  -Y3(0 


A8IYY  =  |Y2-Y3(i)| 
AB2YY=|Y2-Y3(2)| 


"K^ 


A6IYY-A&2YY 


=  0 


1  =  2. 
X3(l)=X| 


<0 


JL 


ERROR 
STATEMENT 


NO  ■  NO  +  I 


i=i 
X3U)=XI 
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A»  1  +  {1Z?/W 


X30)  = 


X3(2)= 


-S-*-'VBa-'«AC 


2A 


.6-V6*^aT 


"ZA 
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A&SXXl=|)(l-X3(l)| 
ABSXXi-|X2.-X3(*)| 


Y3l=V"vl-(X3(l))1 
Y32=  -V31 


-^-(  ABSXXI-A.BSXX2. 


AaSYYt*|YZ-Y3l) 
AfiSVYfc=|  YS.-Y32) 


vl/     =    O 


DH 


ERROR 
STATEMENT 


Y33  =  fv^-(x3W)2- 
y3f=-Y33 


=0  Z'  >>o 

A6SYYI  -A6SYY2. 


ERROR 
STATEMENT 


WO=NO+l 


© 


<0 


JZ 


A<JSYYl  =  |Yi-Y3i| 
AfeSYYZ  =  |Y2.-Y3V| 


^AfeSYYl  -  AfiSYYaj-^t 


JU. 


1-  I 

Ys(r)=Y32. 


*o 


I--2. 
Y3(I)'Y3^ 


ERROR 
STATEMENT 


\/_ 


Y3(J>Y3l 


No*  ft/o +  1 
M*M-H 


I«2_ 
Y3(l)  =  Y33 


© 
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reuspo  ^Ki^-xi^  +  CVaCiVYO1 


RELATIVE 
SPE6.D 


TD  =  m  /RELSPD 

RT  =  PTOT  +■  TOT 
+  CLDT  4-RELT 


XX=0 
SUBT=  RAX/U3 


C    SUBT-RT      J 1 


-» 


>0 


e 


rx=i  , 


REACT  IOW  TIME  TOO 
GREAT.  SU6  CAW 
kTLU  VEHICLE.. 


Tli=TO+RT 
OIST  =  U3  X  TB 


X*/sMST"SINfTHET^ 

+  X2. 
V^OlSTxCOS^KETA) 


\f 


MS0«Ts^W)V(W)* 


JSjl 


8LKVD  TIME 
OISOAT 
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-0 


SIGNAV^PKCWTW  «  RAI 

SlGSRsWNTSxRAI 
Sl6SBft=RAIxSW(e«GER) 


JZ 


SISAPR=Ysi<3SRxSlQSBK 
1      -4-SIGAPR* 


SIG3T 


il£- 


SAVE  =.TB 


c 


J^_ 


TS-TI 


£0 


K^I> 


>0 


ALERTED  CASE 


ERROR 
STATEMENT 


N0«W0+| 


^ 


>0 


ERROR 

swEMorr 


±. 


IFLAG  - 1 
T6*Tt 


W0  =  N0+1 
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\JN=  MO  +-I 


500 
N=l,NuM 


DELTAL'SW 


DELTAL-  OEV.TAL-.3fcO* 


Df  LTAL  »  <*0° 


STft  X  DFLV 
DFlTA(N)  =.  DELTAL 


EVASION 
TURN 


1 


ilL 


DN=  |De\.TA\_-3fcO'| 


£o 


DN=  .ooo<i 


>-*o 


>0 


60 


<0 


DEV-TAu  -3(,0 


5    <0 


dm'Ideltau-ro0! 


DM-.ooo^ 


5 


>o 


DELTAL.-TO 


Ab=0El.TA(N)s^0O 
X(N)  *  STR 

Y(n)=str-h(tb- 

ttt)x  U3 


evasion 

TUR.N 


r 


>Q 


DELrA(N)=:0 
TTTsO 
A6r6BsO 
STRsq 


A6-68^-(lto-oev.TA^ 
DELTA  (*/)  =  A6 
TTT-  STR,xA6 

ITT-  U3 


EVASION 
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X (N )  -  STR.  X  SIN  (C,e>)  +  U3 «  (T6  -TTT)  x  COS  (AG) 

Y(N)  s  STR-  STB.  x  COS  (GG.)  +  U3  x  (TB-TTT)  X  SIKJ  (*&) 


W0-NO  +  | 


XXX  ■  X(N) 
YYY=  YW) 


3Z 


T&  =  SAVE.-TI 


X(N)  «TD  «  U2.  x  COS  (Afi) 
V(^)=16xu2xSWCAB) 


X(«)*XXX  +*(N) 
Y(h)»W*Y(*0 


D6LTAL-OeLTAu  +-  ETA 


XX  =  01  ST 
YY^O 


BIAS  («)  «Y (X(w)-DIST)x-h   (Y(V))' 


I 


D=    IiXA5(N)/siSJT 
AKTLL  =  £XPaN6/SI<ajr 


COMPUTE  PPK(N) 
FOR  ALERTED  CASE 


Compute  evasion    P£o<3acility    op    success  y  phrj 

FOR   TK£   (T1&    CVASXOlJ     usrwo-    THE      ClACOLAQ. 
tOV(T«/V*e      FUNCTION/. 


ate. 


PPk(N)-PK 

SUM  =  SUA1  +  PPK^N) 


iki 


f«oeA6ri.iES 
of  success 


500 


© 
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WTPAV&«6UW\/viM 


A. 


?MG-  Sum/uN 


AVERAGE   EVASION 
SUCCESS    PKOB. 


SVJMUIVU  =  O 


S"0£ 


WEIGHTED     EVASION 

success   fftoft. 


SUMUWU=SUMUWU  +UNU(W) 


>0 


/      SUttUNU  \ 


=  o 


STATEMENT 


WEIGHT  ED    PP>06. 
NOT   USED* 
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W=  (ui  -  AA  (j)  x  UERR)  XTB 


JZ 


GOO 
J- 1,5 


UNALERTED  CASE 


T 


THE.TER=  AA(T)  x  THDEL 


T 


SUBX(I,J) = W  x  cos (thete r) 


-il£- 


K  AS(1 ,  J)  ^fxX-SUBXCljJV+CYV-  SUBY(l^a 


I 


D=  BIf>,S(I,T)/siGJT 
AKILL  =  EXPRNG/SIGJT 


COMPUTE    UNALERTED    PROBABILITY  OF  SUCCESS  ,  PK,FOR,  (I,J)™ 
CELL     OS1WG     CIRCULAR     CO VE RAGE.   FUNCTION  . 


2>L<_ 


Sum  «  Sum  +  PK 


PROG>  =  Sum  x.O* 


-&- 


NO=NO+  I 
M=M+l 


UNALERTED 
SUCCESS     PROS. 
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SUBROUTINE   INPUT 


(    RETURN     k 


SUM=  O 

SUMM=0 

IFLAG=0 


<^D 
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READ  2.G  INPUT 
PARAMETERS 


I 


I 


4 

26 


unu(n),n«i,nun 


f  RETURN    J 
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